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Abstract. We are motivated by the problem of control for a non- homogeneous elastic string with 
memory. We reduce the problem of controllability to a non-standard moment problem. The solution 
of the latter problem is based on an auxiliary Riesz basis property result for a family of functions 
quadratically close to the nonharmonic exponentials. This result requires the detailed analysis of an 
integro-differential equation and is of interest in itself for Function Theory. Controllability of the 
string implies observability of a dual system. 
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1. Introduction. The conditions of controllability for the broad class of the 
linear oscillating structures have been under consideration since the classical papers 
by H.O. Fattorini and D.L. Russell (see [3^, [30], the survey [37]: and book [B] for the 
history of the subject and extended list of references). The method of as well 
as many other subsequent papers, is based on the properties of exponential families 
(usually in the space L^(0,T)), the most important of which for Control Theory are 
minimality and the Riesz basis property. Recent investigations into new classes of 
distributed systems such as hybrid systems and damped systems, as well as problems 
of simultaneous control, have raised a number of new and difficult problems in the 
theory of exponential families (see [I]-[II], [H] , [U] , [2S] ) . 

We say that a mechanical system, e.g. a string, is controllable if, for any initial 
data by suitable manipulation of the exterior forces, the system goes to the given 
regime. According to the classical scheme outlines above, the solution of the control- 
lability problem is based on an auxiliary basis property result. The last result is of 
independent interest and represents the central part of the paper. 

In this paper, we study controllability of an oscillating string the material of which 
has the memory. Such a string is described by the equation of the form 



(1.1) vu{x, t) = y,,(x, t) + N{t- T)y,,{x, T)dT, {x, t) e (0, /) x (0, T), 



and we suppose that a control force is applied to one end of the string. Here and 
below, N{-) is the memory kernel. 

In the short review of the literature below, we discuss the results on both heat 
and wave equation with memory since, formally speaking, controllability for the heat 
equation may be studied in a way similar to one for the wave equation. Also, the 
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results on control for heat equations with memory motivated our study. We consider 
only linear models. Also we slightly change the authors' notations to make them 
uniform. 

To the best of our knowledge, G. Leugering has been the first author to study 
controllability for viscoelastic systems (see [27l[28]). The model studied in [27] leads 
to the initial boundary value problem for the heat equation with a memory term. 
Since the case of the constant coefficients is considered, Laplace transform is used. 
Monotonicity of the memory kernel is essential for the proof of the exact controllability 
of the system in finite time. The existence of the time-optimal control is shown. An 
abstract differential equation of the second order with a memory term is considered 
in and the generator of the appropriate semi-group is studied. 

It has been known for a long time that the infinite speed of propagation exhibited 
by the heat equation is not quite physical. The corresponding modification of the 
Fourier law results in the heat equation with memory (see [H], [M]). 

Controllability for the heat equation with memory is studied in |15l in the multi- 
dimensional case 

(1.2) yt{x, - 7 A y{x, t) - j N{t - s) A y{x, s) ds = f{x, t)xc.{x) , 

Jo 



{x, t) enx (0, T), y{x, 0) ^ya{x), x e n; y{x, t) = 0, (x, t) e 91) x (0, T) . 

Here the symbol denotes the characteristic function of the open set w C fi. The 
controllability is discussed under the strong assumptions on the memory kernel, which 
are alien to the mild assumptions we make in this paper. Specifically, N{-) is supposed 
to be completely monotone, i.e. (— !)•' A^^^^ (i) > Vi > 0, j = 0, 1, .... Under addi- 
tional assumptions on its Laplace transform, approximate controllability is proved 
for T > if 7 > 0. The one-dimensional equation on x G (0, Z) is studied under the 
similar assumptions for the memory kernel. The exact controllability is proved for 
T > l/yjN{{)) if 7 = 0. 

In |41| , a multi-dimensional problem with more general memory kernel is consid- 
ered 

(1.3) yt{x,t)-\j-i N{t- s,x)\j y{x,s)ds ^ f{x,t)xuj{x) , 

Jo 

subject to the similar initial and boundary conditions. The memory kernel at i = s is 
strictly positive (more precisely, mg < A(0,a:) < mi, Vx G fi.) Also some geometric 
conditions on the set uj are assumed. Then there exists Tq > such that the system 
is exactly controllable for all T > Tq by means of control / € L^[(0,r) x w]. The 
method is based on Carleman estimates. 

The paper [42] represents a generalization of |41j for anisotropic and nonhomoge- 
neous medium, i.e. the equation 

(1.4) yt{x,t)-^ (a'^ix) J A^(i - s, x)y,. (x, s) ds") ^ fix,t)xUx) 

is considered. The similar results on exact controllability and observability are proved 
under the geometric conditions, which are more complex than in [41j. On the other 
hand, if the geometric conditions are not satisfied, a negative controllability and 
observability result is presented. 
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If we consider (|1.2p with 7 = 0, then differentiating (ll.2p with respect to time for- 
mally yields the wave equation with memory. Yet, generally speaking, controUabihty 
of the heat equation does not imply controllability of the wave equation. 

The model 

(1.5) ytt=NQAy+ f N{t ^ s) A y{x, s) ds - ayt 

JQ 

is considered in [3T]. Here A^o > 0, < 0, and a > 0. The exponential decay of the 
solution is proved assuming the memory kernel N{t) ^ and decays exponentially. 
The reachability problem for the wave equation with memory 

(1.6) Utt -Uxx+ / e"''(*~''^Ma:a:(a;, s) ds = 

^0 



is studied in |31) . We emphasize the special type of the memory kernel in (jl.6p . 
The goal of the current paper is to prove controllability of a similar system with an 
arbitrary memory kernel (except smoothness requirement). 

In the next group of papers, the memory of the material is not necessarily taken 
into consideration. Yet, these papers contain some technical ideas which are important 
to us. The wave equation with the time dependent tension is studied in [1], [2], [4], 
and [5] . The proof of Riesz basis property of a family of the time-dependent functions 
in [T] and [3] is based on the assumption that the tension varies "slowly enough" with 
time. This restriction is removed in [5j. To the best of our knowledge, the papers 
[I]; [S] represent the first attempt to apply the method of moments to equations 
with time dependent coefficients. The new difficulty in this case is the absence of an 
explicit representation for a family of functions arising in the moment problem. This 
fact, which substantially complicates the analysis of controllability, is common to the 
heat equations with memory. 

The papers [Ml |3S] are especially important to us. 

The one-dimensional version of (|1.2p with 7 = and without any monotonicity 
assumptions on N is considered there: 

(1.7) yt{x,t) ^ ay{x,t) + N{t - T)yxr,{x,T)dT . 

Controllability problem is reduced to a moment problem with respect to a family of 
functions. Techniques developed by L. Pandolfi in ^34l 135) allow proving that the 
aforementioned family forms a Riesz basis in a proper space, and that allows 
solving the exact controllability problem. In |12[ I13j the problem of simultaneous 
temperature and flux controllability for heat equations with memory was studied. 

Similarly to [Ml |3S1 [121 US], we reduce the study of controllability to investi- 
gation of the Riesz basis property of the auxiliary functions, which we call "quasi- 
exponentials" . However, there is a serious difference between the models discussed in 
these papers and our model, (a) In these papers, the transformation 

(1.8) y^e^^y with 61 = -iV'(0)/A^(0) 

leads to an equation of the same form as (1.7) but with A^'(O) — 0. The last condition 
appears to be very helpful for the technical purposes. For the wave equation with 
memory, the similar simplifying condition has the form A^(0) = which can not be 
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achieved by a change of variables preserving the structure of equation. This fact 
sufficiently complicates the analytic problems we solve in the current paper. That is 
where the difference between our model and the models considered in |M1 [351 HH US] 
is essential, (b) We note also that the different techniques we use here results in the 
less restrictive requirements on the kernel N{t) (see Assumption 2 below). 

The present paper is the first work where the boundary controllability of a general 
Id-wave equation with memory is studied. Sharp controllability time is established 
under mild regularity conditions for the (variable) coefficients of the equation. 

The inclusion of the memory term in the model makes it necessary to analyze a 
non-standard moment problem associated with the solutions to a family of the Cauchy 
problems for an integro-differential equation. We use diverse methods of Functional 
Analysis and Asymptotic Analysis to prove that these solutions (we call them quasi- 
exponential functions) possess the Riesz basis property in a proper space. This 
property allows us solving the moment problem and is of interest in itself for Function 
Theory. 

The paper is organized as follows. In Section 2 we discuss the statement of the 
problem and the auxiliary Sturm-Liouville problem. In Section 3 we derive a series 
representation for the solution and the moment problem with respect to the quasi- 
exponential functions. Section 4 is the central in the paper. Here, we prove basis 
properties results for these quasi-exponentials. In Section 5 we solve the moment 
problem and prove the main theorem about controllability, i.e. give the conditions of 
exact controllability of the string with memory. We also prove the observability of 
the dual system. 

2. Statement of the initial boundary value problem and controllability 
problem. The auxiliary Sturm Liouville problem. For any T > 0, we consider 
the following initial boundary value problem (IBVP) for the hyperbolic type PDF 
with the convolution term 

(2.1) p{x)ytt{x,t) = {Ay){x,t)+ f N{t ~ T){Ay){x,T)dT, {x,t) e (0,1) x {0,T), 

Jo 

(2.2) y,{0,t)^f{t),y{l,t)^0, 

(2.3) yix,0)=yt{x,0)=0. 
Here A is the differential expression 

(2.4) iAip){x) = iaix)ip'{x)y + /3(x)^(a:) . 

Above, p{x) is the density of the string; I its length; a{x) its modulus of elasticity; 
N{-) its memory kernel; the function f{t) is control, i.e. the exterior force acting on 
the left end of the string. 

The IBVP (I2.ip ~ (l2.3l) describes the small (linear) transverse oscillations of the 
string with memory. Since the similar IBVP without memory term has a unique 
solution if the functions p, a, b, and/ are smooth enough (j26j). we may hope that 
the integral term in will not change this fact if the kernel N is smooth enough. 

This paper solves the exact controllability problem for the system (|2.1[) - (|2.2|) . 
More precisely, for any T > 0, we describe the reachability set 

|(/(.,r),y/(.,r)): /eL^(o,r)| 
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of the system. 

We proceed further under the following 
Assumptions. 

1. p and a are strictly positive functions from C'^[0,l], (3 e C[0, . 

2. N" e L2(0,T). 

Remark 1. Assumption 1 is used to claim the well-known results on the aux- 
iliary Sturm-Liouville problem (see Lemma 1 below). Using recent results on the 
Sturm-Liouville theory (see, e.g. )40|). we may weaken the regularity assumptions: 
p", a", b € L\0,l). 

We further introduce the functional spaces, 

(2.5) H = LliO,l), Hi = {^eH\0,l),^(l)^0,} 

where H^(0,l) is the standard Sobolev space. Then there is a natural continuous 
embedding Hi C H, which leads to the natural embedding of H into the dual space 
v.-! := (Hi)' to Til. We thus obtain the triple of spaces (so called rigged space): 

Hi c 7^ c T~L—i ■ 

We now may define the operator A as the differential expression (|2.4I) on the space T-L 
with the domain H'^(Q, l)n'Hi. 

We consider the Sturm-Liouville problem associated with (I2.4p . (|2.2p : 

, , (a{x)(p'{x)y + f3{x)ip{x) + Xp{x)ip{x) =0, < x < / ; 

^ ' ^'(o) = ^(0 = o. 

The following result is well-known (see [E], Ch. X; [43], Ch. VI; [29]; [32]; [33]; 
[22], and the references cited therein). 

Proposition 1. Let Assumption 1 be satisfied. Then the Sturm-Liouville prob- 
lem i2.6\) has the following properties. 

(a) Its spectrum {Xn}'^=i is pure discrete, simple, and real with the only point of 
accumulation at +oo. 

(b) The asymptotic representation holds, 

(2.7) Wn := a/A^ = a/A5[+ as n -> oo . 

Here 



(2.8) ^,.l(n-l).n>l^. L^jje^,.. 

(L is the optical length of the string, i.e. the time of propagation of the waves along 
the string). The eigenfrequencies w„ are separated, 

(2.9) inf \ujn-iOk\>0. 

(c) The corresponding eigenfunctions {^n}^=i form an orthogonal basis in % 
(which we assume to be orthonormal) , 

(2.10) / p(x)tpn{x)tpmix) dx = 61^ , [Alfn, (fm) L^OX) = S'^ . 
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(d ) The estimates hold 

(2.11) |<^„(0)|xl. 
This relation means that 

< inf \(pn{0)\ < sup \(pn{0)\ < oo . 

new „gN 

Proposition lb imphes 

Proposition 2 _6:, Sec. II. 4], [39]. Put Tq = 2L. The following statements are 
valid. 

(a) If Xn 7^ for all n, i.e. a;„ ^ for all n, then the family £ = {e^''^"*}„eN 
forms a Riesz basis in L^(0, Tq). 

(b) If Xn =0 for some n, then the statement (a) remains valid for the family 

f = {e±^""*}„:<.„#oU{l}U{^}• 
('cJ For T > Tq, there exists an inhnite family of exponentials £j such that the 

family £ [J £f forms a Riesz basis in L^{0, T) . 

(d) For T < Tq, there exists an inhnite subfamily £^ C £ such that the family 

£ \ £T forms a Riesz basis in L^{0, T) . 

All the statements of Proposition 2 remain valid if we replace the family £ by 

|g±icd„ with any e C, since the operator of multiplication to e^^ is bounded 

and boundedly invertible in L^(0, T) and therefore preserves the Riesz basis property. 
Remark 2. According to Proposition 1 (a), only a Enite number of A„ may be 

negative. In this case, the corresponding aj„ are pure imaginary, ujn = ivn, 3,nd for 

dehniteness, we assume > . 

The following result specifies the regularity of the solution to the IBVP (|2.ip - 

(1231). 

Theorem 1. Let Assumptions 1 and 2 be satisfied and f E L^{0,T). Then 
2/eC([0,r];Hi) andyt€Ci[0,T];n). 

Below we use the notation y-^ for the solution to the IBVP subject to the (con- 
trolling) force /. 

The following results states the controllability property of the IBVP (|2.1[) - (|2.3|) . 
Theorem 2. Let Assumptions 1 and 2 be satisfied. Then 

(a) For T > Tq, the system i2.1\) - ^2.3\) is exactly controllable, i.e. 

(2.12) |(y/(.,r),y/(.,r)) : / e ^^(0, T)| = Hi x H . 

(b) For T < Tq, the system ^2.1\) - ^2.3\) is not approximately controllable, i.e. 

(2.13) cl |(y/(-,r),y/(.,r)) : f E LH0,T)'^ ^ Hi x H . 
Moreover, 

(2.14) codimcl |(/(-,T), y/(-,T)) : /e L2(0,T)| in Hi x H = (X) . 

Remark 3. The results of Theorems 1 and 2 can be easily extended to other 
boundary conditions and other types of boundary control and source-type controls 
like g{x)f{t) with the given g. 
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When this paper was ready to submission, we have been informed that the result 
close, though less general, to Theorem 2(a) (controllability of the system described 
by equation (|l.ip for T > 21) was obtained in the paper by P. Loreti, L. Pandolfi 
and D. Sforza, "Boundary controllability and observability of a viscoelastic string" 
(submitted). 

3. Moment problem. As we show below, the moment problem originates from 
the terminal conditions 

(3.1) y{x,T) = vo{x), yt{x,T) = vi{x), x€[0,l], 
where 

(3.2) vqEHi, vien. 

First, we derive an appropriate representation for y{x,t). For that, we consider the 
inner product (Ay, (y5„)^2(Q and integrate by parts twice to find 

(3.3) / ((ay')' + Py)^ndx = K„/(t) - A„ / py^pn dx . 
Jo Jo 

Here we use the relations (|2.2p and (12.61) and introduce the notation, k„ = — a(0)(/3„(0). 
According to (pHI) . |k„| x 1 . 

We further expand the solution y{x, t) to the IB VP (|2.ip - (|2.3l) into the series with 
respect to the basis of the eigenfunctions (pn (see Proposition 1) 

(3.4) y{x,t) = ^ an{t)Lpn{x) . 

n>l 

We then multiply both sides of the basic Eq-n (j2.ip by an arbitrary eigenfunction (fn , 
integrate over (0, Z), and use the identity (|3.3p . We finally find the integro-differential 
equation for the coefficients a„(here and below we revert to the notation a;„ = \/\^), 

,2 „ l4-\ i , ,2 



(3.5) a„(t) +a;>„(t) +a;^ / 7V(t - T)a„(r)dr = k„ g(i) , 

Jo 

subject to the homogeneous initial conditions (see (|2.3|) ') 

(3.6) a„(0) = 0, a„(0)=0. 
Here we introduce the new "control" g(t), 

(3.7) g{t) = f{t)+ f N{t-T)f{T)dT. 







We notice that the map f ^ g \s linear bounded and boundedly invertible in L^(0, T). 
The following estimate and equality are standard (see, e.g. (6] Sec. HI.l]) 

(3.8) ll2/(-,i)ll?., ll2/*(^0llw = E 

n>l n>l 

Our next goal is to obtain a representation for a„(i). For that, we introduce two 
families {c„(t)}, {s„(t)}, n > 1 of solutions to the homogeneous integro-differential 
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equation corresponding to (I3.5p . satisfying different initial conditions, 

Cnit) + ojI c,,{t) + ojI f N{t ~ r)c„ (r)dT = 0, c„(0) = 1, c„(0) = , 

(3.9) ^0 1 n > 1 

Sn{t) + S„{t) +UjI N{t ~ T)s,,{T)dT = 0, S„(0) = 0, S„(0) = LO,, . 

Jo 

(In the absence of the memory, N{t) = 0, it follows c„(i) = cosa;„t, s„{t) = sincLi„t.) 
Then the solution of the Cauchy problem p.9|) has the form 

(3.10) a„{t) = K„ f 9{r f''^^^^^ dr. 

Jo 

From p.9p it is easy to derive that s„(^)/a;„ ~ Cn{t). Therefore, 

(3.11) d„(t) = K„ / g{T)cn{t - t) dr . 

Jo 

Hence, we have found a representation p.4p . p.lOp . p. lip for the solution y{x, t) and 
its time derivative yt{x,t). Then the Eq-ns (j3.1|) can be written in the form 

(3.12) ^ an{T)ipn{x) = va{x), ^ d„(r)(p„(x) = vi{x) . 

n>l n>l 

We further use the representation p.lO[) . p. lip and the orthogonality condition (|2.10p 
to conclude 

g{T) Sn{T -T)dT = W„Wo.n, 

(3.13) -^0^ ri>l 

f^n / 5(r)c„(r - r) dr = t;i,„, 
Jo 

where we let 

(3.14) Vo^n = (Vo, 'Pn)n, Vl.n = {vi,^Pn)n ■ 

Inclusions (13. 2p imply that {u}nVo,n} , {'^i.ti} G 

The equations p.l3p represent the desirable moment problem: Find g G L^(0,T) 
given {uo,n} , {vi^„} ■ 

Theorem 1.2.1 of jGj implies that the solvability of the moment problem p.l3p is 
the direct consequence of the basis property of the system {e^(t)}. First, we discuss 
the basis property (Section 4) and then solve the moment problem (Section 5). 

Instead of the family {c„} U {s„} we may consider the family £ ~ {e^(i)}i°, 

e^{t) = Cn{t) ± iSn{t) 

that satisfies the equations 

(3.15) et(t)+^let{t)+ujl [ Nit - T)et{T)dT ^ . 

Jo 

The initial conditions for the functions (t) originate from the corresponding condi- 
tions for Cn{t) and Sn{t), 

(3.16) e±(0) = l, e±(0)=±*wn. 
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If for some n, uj^ = 0, there is only one solution to the Cauchy problem (|3.15l) - (|3.16p . 
e„ — 1. In order not to "lose" one solution, we proceed as follows. We note that 
the Eq-n (I3.13P of the moment problem contains the ratio, s„(-)/aj„. The expression 
Sn{t)/ujn should be understood in the limiting sense if w„ = 0. We show that for any 
finite 

(3.17) lim '-^=t. 

Indeed, the initial value problem p.9|) for Sn{t) may be easily transformed into the 
following integral equation 

(3.18) s„ (i) = sin cj„t — cj„ / sinaj„(t — r) / N{t — z)sn{z) dz dr . 

Jo Jo 

The limiting representation p.l7|) may be now easily derived from p.l8|) . Therefore, 
if ujn — for some n, then we put the corresponding e^{t) in the family {e^(t)} to 
be e+ = 1 and e:^ — t (see Proposition 2b). 

4. Basis property o{ £ — {e^}. In this section we develop new technical ideas 
that generalized the methods used in our previous papers [U |4] and in papers by L. 
Pandolfi [MlES] (see also our joint paper [5]). 

We will use two definitions. 

Definition 1. A sequence {e„} is said to be uj— independent in a Hilbert space 
H if tlie conditions 

{C„} € and = 

imply that ^„ = for every n. The convergence of the series is understood in the 
norm of H. 

Definition 2. A family {e„} is said to be an C-basis in a Hilbert space H if it 
forms a Riesz basis in the closure of its linear span. 

The following theorem represents the main analytic result of this paper. 
Theorem 3. Let Assumptions 1 and 2 be satisfied. Then 

(a) the family £ — {e^{t)} forms a Riesz basis in L^{0, Tq) ; 

(b) for T > To, there exists an infinite sequence {ilin} of real numbers and the 
corresponding family f J of solutions to the Cauchy problems i3.15\) - f3.16\) such that 
the family £ [J £j forms a Riesz basis in L^{0, T) ; 

(c) for T < Tq, there exists an infinite subfamily £^ C £ such that the family 
£ \ £T forms a Riesz basis in L^{0,T) . 

Our proof is based on the Bari theorem [W, f23', p. 317], O Theorem 15 p. 38]. 

Theorem (Bari). Any uj— linearly independent sequence {e^} which is quadrat- 
ically close to some Riesz basis {e„} in a Hilbert space H, is itself a Riesz basis. 

From Proposition 2 it follows that the family {e^*"^" t+Af(o)t/2| fgrnis a Riesz basis 
in L2(o,2To) and an /I-basis in L^{0,T) for T > 2Tq. Therefore, two steps suffice to 
prove the basis property of the family £ = {e^(i)}. 

1) We find the asymptotic representation of the family £ — {e^} to prove that it 
is quadratically close to the family {e±»"" *+^(o)*/2} in L'^(Q,T). 

2) We prove that the functions e^{t) are w— linear independent in L?{Q,T) for 
T> To. 
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4.1. Quadratic closeness. At the first step, we may assume t ^ [0,T] with an 
arbitrary chosen T > 0. We construct the leading term of the asymptotic representa- 
tion for e^{t) as n — >■ oo. That allows us to prove the quadratic closeness. Our proofs 
are mostly based on numerous and cumbersome estimates with the main technical 
tool to be the Gronwall inequality along with the direct asymptotic estimates. 

The integro-differential Eq-n p.lSp subject to the initial conditions (|3.16p is 
equivalent to the integral equation 



(4.1) e±(t) = e^^*"^"* -a;„ / sina;„(t-r)/ N{t - s)e^{s) ds dr 

Jo Jo 



= e^'""'-a;„ / sina;„z / N{t ~ z - s)e^{s) ds dz 
Jo Jo 

or after integration by parts 

(4.2) e^(i) = e=^*""* - /" N{t - s)e^{s) ds + N{0) [ cos cUnZ e^{t - z)dz 

Jo Jo 



t rt-z 



+ / cosa;„z / N' {t — z — s)e^{s) ds dz . 
Jo Jo 

We rewrite the last integral equation as follows, 

(4.3) e± (t) = e±'""* + /* K [t - s)e± (s) ds 

Jo 

where the following kernel N*{-) is introduced: 

(4.4) N*(t) = -N{t) + N{0)cosu}nt+ [ cos uj„z N' {t - z)dz . 







Below, we will need to analyze the asymptotic representation to the solution of a 
slightly more general integral equation than (j4.3p 



(4.5) (i, /i) = e±-"* + / e-'^^*-^) (t - s) [s, a.) ds 



with some parameter fi. Obviously e^(t) — v^{t,0). 

Obviously the kernel e'"* N*{t) is bounded provided N' £ L'^iO.T), 

(4.6) max \e-^'* N*{t)\ < ((l + ci)max \N(t)\+ci max / |iV'(z)|dz)el^l'^ = ca, 

[0,T] ^ [0,T] [0,T] Jo 

ci = max I cosa;„t|, t G [0, T] . 

n 

(We remind here that some first w„ may be complex). More than that, N* e C[0, T] . 
In (|4.6p and below, c with a lower index denotes a positive constant that is independent 
of both n and t. 
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Eq-n (j4.5p and the Gronwall inequality now imply 

(4.7) \vtit)\<C3, te[0,T]. 

This result, along with Eq-n (14.51) . implies that v^{t) e C^(0,T); for the function 
e^*'^"*, this is also seen from p.lSp . 

To estimate v^{t) more precisely, we rewrite (|4.5I) as follows 

(4.8) vtit) = + (AvtW) + [B^vtm + {C,,vtm 
where we introduce three integral operators 



Jo 

{Br,vt)(t) = N{Q) I e~^'^'-'UosuJn{t-s)vt{s)dz- 







{Cnvt){t)= cosujnz 3- N' {t - z - s)v^ {s)ds dz . 

Jo Jo 

We proceed with the estimate of v^{t) in two steps. 
Step 1. The formal asymptotic representation for large n. 
First, we find 

Jo Jo 



* sinujnZ 







dz . 



iV'(0)e''^ v^{t -z)dz+ / e-^(*-^) N"{t - z - s)vt{s) 



so that 

(4.9) |(C„^;±)(^)|<^ if N"eL\0,T). 

n 

Everywhere below, the symbol O means an estimate that is uniform with respect 
to t e [0,T]. 

Since we proceed formally at this moment and are interested in the leading term 
of f^(t) as n — 00 only, we remove the term (C„ v^){t) = O(^) from the analysis, 
hence considering the (approximate) integral equation 

(4.10) = e±-"* + {Avtm + (B„ vtm . 

The solution of (I4.10p may be constructed by the (convergent) iteration series 

(4.11) vt{t) = {I + A + A^ + ... + Bn + Bl + ABn + B^A + ...)e±*""* . 
We have \A^ e=t''^"*| < C5/71. Indeed, 

(4.12) (Ae±^""*)(i) = - / e-"^'-'^ N(t~ s)- . 
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Integrating by parts yields the desired estimate provided N' e L^{0, T) (Assumption 
2 is even stronger, N" e L^(0,T).) The estimate for the higher powers of A is then 
obvious. We further evaluate the term (B„ e^'""*)(t) 

(B„ e±''^"*)(f) = iV(0) / e-^(*~'') cosa;„(i - s)e^"^-' ds = 







2/i \n 
Repeating the previous calculations yields 

(4.14) (Bl e±*""*)(i) = e±*""*-''* (e''* - 1 - m*) + O (^^ 

According to (|4TT2)) and (I4TT3)) . the term (ylB„ e=^*'^"*)(t) is of the order 0(i) (we 
again assume N' £ L'^{0,T).) 

We conclude that the only terms in the iteration series (|4.1ip which produce the 
leading terms for v^{t) are generated by the leading terms of (|4.13p . (|4.14p . etc. 

(4.15) v^{t)^{I + B^ + Bl + ...)e^^--*. 

The leading term of the asymptotic representation for (i?^ e^'""*)(f), k > 2 may 
be found with the help of the consecutive applications of the operator Bn and removing 
the terms of the order O(^). We finally get 

(4.16) v^it)^Y.(B'ne^'^'''m = ^^'^^'-''i: E ^-^(^Y 

j>0 M>0 j>M ^ ^ ^ 



E 



_ ±tu>„t-tit ^ I ^(0) „N{0)t/2 _ at/2 \ , 

(iV(0)/2^) - 1 1^ 2m ' J + - 

We need two specific cases of this formula. 

Case 1. If — >■ 0, the integral equation (|4.5p becomes the equation (|4.3p . Hence, 
the representation for ej(0)t has the form 

(4.17) e±(t) = hm v^{t,^l) ^ e±-"*+''* + ... 

where we introduce the notation which appears below on the regular basis and will 
be used in the exponents only 

(4.18) i/ = 7V(0)/2. 

Case 2. If /i A^(0)/2 = ly, the integral equation (|4.3I) has the form 

(4.19) z;±(t,;.)=e±'""*+ f e'''^'-^^ Kit ~ s) vt{s.v) ds 

Jo 
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and the following representation holds 

(4.20) vtit)^ lim t;±(t,M) =e±'""*(l + 7V(0)</2) + ... 

The correction terms in (|4.17l) and (|4.20[) seem to be of the order 0(i) though 
that has not been proved yet. 

We now briefly consider a particular case of the memory kernel N{t — s) = 
Qg-'K'-s) (see [31]) where the functions e^(t) may be found explicitly (and hence 
asymptotically). Indeed, if we apply Laplace transform to the problem 

Jo 



e±(0) = l, e±(0)=±ia;„, 

we find 

" ^ ^ 2^iJc P^ + VP"" + l^Ip + u:Ua + v) 

with the usual contour of integration C in the inverse Laplace transform. An elemen- 
tary perturbation theory allows finding the poles of the integrand 

, ■ 

pi,2 = ±iU}n + 2 + P3 ^ -a. - V + ■■■ 
Thus, the asymptotic representation for e^{t) is the same as in ()4.16p 

e±(i) =e±^'""*+''* + of — 

where we use the fact that N{0) = a. 

Step 2. Justification of tie asymptotic representation ^4.17]) . 

Lemma 1. Let Assumptions 1 and 2 be satisfied. Then tiie family £ = {e^(t)} 
is quadratically close in L^{0,T) to tlie family {e='=*""*+'^*}. 

Proof. We actually need to prove that the order of the correction term in (|4.17p 
is 0(i) uniformly on t G [0,T]. For that, we introduce the new function, 

(4.21) i?±(t)^e±(t)-e±-"*+''*. 

First of aU, E^{t) G C^[0,T]. We now derive an integral equation for E^{t). For 
that, we substitute e^(t) = e=^'""*+''* + E^{t) into Eq-n (|431) . After some long but 
elementary calculations, we find 



(4.22) E^{t)= N*it-s)E^{s)ds + F^{t), 

Jo 

where we let 

(4.23) F^{t) = - f N{t - s)e±''""^+'^^ ds 

Jo 
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^(0) 



7V(0) ± Aiuj„ 



+ [ COSUJnit-T) [ N'{T-s)e^"^"'+''-'dsdT. 

Jo Jo 



Integrating by parts in the underlined integrals yields 

^(0) „±..„t+.t ^ Njt) 



(4.24) F,f{t) = -—-^^-—e 



±iujn + N{0)/2 ±iuj„ + N{0)/2 

* N'(t - S) e'^'""''^''" ^ , ( ±^u^^t+^t _ Ti'^r.t^ 

^ ' ±iLOn + A^(0)/2 ^ 7V(0) ± 4iw„ r J 



* sin cjji — r) 











so that 



(4.25) F±(t) = ^ where |0±(i)| < if N" e L^{Q,T) . 

n 

Applying the Gronwall inequality to Eq-n (|4.22p and using the estimates for the kernel 
(j4.6p and for the free term (|4.25p we find 



(4.26) max |i?±(OI<-- 

te[o,T] n 

Hence, according to (|4.2ip . the system {e^{t)} is quadratically close to the basis 

The proof of w— linear independence of {e^ (t) } below requires analyzing the prop- 
erties of {Fjf{t)y and (F,^(i))". Representation (I4.24p for F^{t) implies the uniform 
with respect to i G [0, T] asymptotic representation for its derivative 

(4.27) (F^it))' = - Ar(0)e±'""*+''* + L±^'^,^t+ut ^ ^^^u^t 



+ / cosw„(t-r) / N"{t- s)e^''^^'+'''dsdT 
Jo Jo 



N'{Q) 
' NiO) 



(e^^* - 1) +i/;i,„(i), where Vi,n(t) =o{ — 



Here we explicitly show all terms of the order 0(1) and introduce the correction 
term, ipi,nit)- Representation (4.18) shows that if a;„ = for some n, the corre- 
sponding is independent of w„. In this case, we let Vi.n — 0(1). The similar 
representation for the second derivative has the form 

(4.28) {F^{t))" ^ -N{0){±tu;^ + ^)e±-"*+'^* 
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-w„7V'(0) / sina;„(t-T) [ N" {t - s)e^'^^''+'"' ds dr 
Jq Jq 



Before proving w— linear independence of {e^{t)} we need some auxiliary results 
on the properties of the functions {t) and (t) . 

Lemma 2. Let Assumptions 1 and 2 be satisfied. Then 

(a) i^±(0) = 0, F±(0eC2[0,T], 

(&) |(F±)'(t)| < Cr, < C8^n, 

^ ^ ^ (c) A^*(0) = 0, £;±(0) = 0, \Et{t)\ <f<c„- 

(d) \{E^nt)\<c,o, \{Et)"{t)\<c,,co^. 



where we let ujn — \uJn\ if iOn ^ 0; uJn — ^ if w„ = . 

Proof. The results in (a) immediately follow from the representation (I4.23[) 
and (|4.28p . The first two results in (c) immediately follow from (|4.4I) and (|4.22p 
correspondingly. The third estimate in (c) is the direct consequence of (|4.26l) . The first 
estimate in (b) is a consequence of (|4.27l) . The second estimate in (b) is a consequence 
of (|4.28p . If ujn = for some n, the corresponding {F^){t) is independent of w„ (see 
(|4.23l) '). To prove (d), we differentiate Eq-n (|4.22p and integrate by parts to find 

(4.30) (E^nt)^ f\N*nt-s)Et{s)ds + N*{0)Et{t) + {F.tnt) 
Jo 



= -N*iO)Et{t)+N:it)Et{0)+ f N:{t-s){Et)'{s)ds + {F^)\t). 

Jo 

Since A^*(0) = -B±(0) = (see (c)) we find the final form of the equation for 

{E^y{t) to be 

(4.31) {Et)'{t) = f K{t-s){Et)'{s)ds + {F^)'{t). 

Jo 

The unique solvability of this Volterra integral equation in C[0,r] follows from the 
continuity of the free term {F^y{t) , see Lemma 2 (a).] Applying the Gronwall in- 
equality to Eq-n (|43T|) and an estimate of \{F,fy{t)\ from (jT29l) (b) yield 

(4.32) |(£;,t)'(i)l<cio. 
Differentiating (|4.3ip yields 



(4.33) (Etnt) = f (KYit - s)iEtns) ds + {F,tnt) 

Jo 
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= ^Kit - s)iEiyis)\lzl + f Kit - s){E^ris) ds + {F,^nt) 



Since N*{0) = (see Lemma 2 (c)) and (£^^)'(0) = (^;^)'(0) (which follows from 
Km ) we find 



(4.34) iE^rit)^N*it){F,fy{0)+ / N:it-s)iE^ris)ds + iF,frit). 

Jo 

The unique solvabiUty of this Voherra integral equation in C[0,T] follows from the 
continuity of the free term {F^)"{t), see Lemma 2 (a).] Applying the Gronwall 
inequality to Eq-n (|4.34l) and using the estimate of \{F^)"{t)\ (see (|4.29l) (b)), we find 

(4.35) \{Et)"{t)\<cnuJ,,. . 
4.2. Lu— linear independence. 

We proceed in two steps, which are similar but not identical to ones in [51135) (see 
also the discussion in Introduction). 

The next result describes the asymptotic structure of the sequence {{E^y{t)}. 
Lemma 2(d) only shows that it is bounded. 

Lemma 3. TJie sequence {{E^y{t)} admits the following asymptotic represen- 
tation 

(4.36) (Etyit) = Di e±»""*+'^* + D2 e±*'^"*+''* • N{0)t/2 + q^^i{t) 

where the sequence {qn,iit)} G uniformly on [0,T] and the constants Di and D2 
are independent of n. 

Proof. The free term of the integral equation (|4.3ip is the linear combination of 
the following terms, 

(4.37) e^"^-\ ^±^^.t+.t^ Vl,n(t), 

and the integral term 

(4.38) In{t)= I cosuJn{t-T) f N" {t ~ s)e^^'^-'+'"' ds dr . 

Jo Jo 

Here the sequence {ipi,n} & uniformly on i g [0,T]. In the integral /n(i), we 
substitute t — s 1-^ s and change the order of integration. After that, the interior 
integral may be evaluated. We find 



(4.39) In{t) - / N"{s)e 



t 



±iuj„t ^ ~ ^ \ o{ — 



ds 



where the last O is uniform on [0, T]. Assumption N" E L^{0, T) implies that /„(i) € 
P uniformly on (0,T). We conclude that the free term is the linear combination of 
three terms, a sequence from P uniformly on (0,T) (see (|4.38p and (|4.39p ) and two 
exponentials (see (|4.37p ). Hence the solution {E^y{t) may be represented as a linear 
combination of the corresponding solutions, which we denote as u^j{t), j = 1,2,3. 
We denote the proofs of the estimates for u„., as Steps El, E2, and E3. 



CONTROLLABILITY OF AN ELASTIC STRING WITH MEMORY 



17 



Step El. First, we show that if the free term of the integral equation (I4.3ip 
belongs to P uniformly on [0,r], then the solution has the same property. Indeed, 
the solution M„,i(t) of the integral equation 

may be estimated as follows 

l<iWI<k™WI + ci / l<i(s)|ds. 

Jo 

According to the Gronwall inequality 

KiWI < \rn{t)\ e^^* < e=^^ \r„{t)\ so that {^^(i)} € l^. 

Step E2. The integral equation with e^*"^"* in the free term coincides with the 
equation (14. 3p . Hence 

Step E3. We use the results of Lemma 1 here. We consider the integral equation 

(4.40) utsit) = + /* N^it - s) u^Js) ds , 

Jo 

rewrite it as 

^sWe-' = e^'-"' f Kit - s) e-'-f*-^) ui,{s)e-'^ ds , 
Jo 

and introduce the new unknown function, v^{t) = ^{t)e~'^^ , which satisfies the 
equation (14. 19^ . Its asymptotic representation for v^{t) is given by Lemma 1, see 
((¥?20| . We conclude 

(4.41) M±3(i) = e±*'^"*+'^* (1 + N{0)t/2) + ... 

The justification of the last asymptotic representation follows one in Lemma 1. We 
introduce the function 

(4.42) Un{t)^ = ui^it) - e±''^"*+''* (l + iV(0)i/2) 
which satisfies the integral equation similar to (|4.3p or (|4.22p 

(4.43) C/±(i) = r Kit - s) C/±(s) ds + G±(t). 

Jo 

For the free term G^{t), after some long but elementary calculations, we find 

(4.44) Gt(.t) = 6^*""*+''* - 6^*""*+''* (l + iV(0)i/2) + 

A^(0) ^ cosunit- 3)6^"^-'+"'' (l + iV(0)s/2)ds + oQ^ =0 
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Repeating Step 1 yields {C/„(t)} G P. This shows that representation (I4.4ip actually 
holds with the error term from and this completes the proof of representation 

Km . . 

The next result describes the asymptotic structure of the sequence {{E^Y'{t)}. 
Lemma 2(d) only shows that it may not grow faster than ii„. 

Lemma 4. The sequence {{E^)" (t)} admits the following asymptotic represen- 
tation 

(4.45) iEt)"{t) = u;„(i^4 e±'"""*+'^* + 6^'""*+'^* + De e^*'""* N{0)t/2 + g„,i(<)) 

where the sequence {qn,2{t)} € l'^ uniformly on [0,T] and the constants are inde- 
pendent of n. 

Proof. The proof is similar to the proof of Lemma 3. We consider the integral 
equation (ITM)) and analyze the free term in it, N*{t){F^)'{Q) + {F^)"{t). The direct 
calculation shows that {F^)'{Gi) = ~N{Q)/2 and also 

(4.46) {F^)"{t) - ±iuJnN{Q)(^- ^ e±*'""*+'^* - ^ e^'""* ± 



smujn{t-T) f N"{t- s)e^''^"''+'"'dsdT^ 



. N'{0) 

As in the proof of Lemma 3, we conclude that the free term in (|4.34p is the Unear 
combination of the following terms, 

(4.47) g±-„t+.t^ gT-„t^ 

where the sequence {ip2,n} G uniformly on t e [0,r]. The remaining part of the 
proof is the same as in Lemma 3. • 

Lemma 5. Let Assumptions in Section 2 be satisfied and T > Tq. Then the 
family £ = {e^(t)} is uj— independent in L^(0, T) . 

Proof. The proof consists of several steps. 

Step A (Linear independence.) We consider an arbitrary finite subfamily 
{e^(i)}„<„f, and demonstrate that the equality 

"0 

(4.48) ^^a±e±(t) = 

± n=l 

is not possible for a nontrivial finite sequence {a„},i<„o. Differentiating the represen- 
tation (j4.ip for e^{t) implies (we note that A„ — uj^) 

(e±(t))' = ±zcj„e±*'^"* + A„ / cosuj„it^T) f N{t ~ s)et{s)ds dr ; 

Jo Ja 

(4.49) (e^it))" = -Ke^"^-' -ul f sinw„(t-T)/' N{t - s)et{s)ds dr 

Jo Jo 



+\n I N{t-T)et{T)dT = -K{e^{t)- I N{t~s)et{s)ds^ 
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Using the last representation and differentiating the identity (I4.48p . we find 

no no 

(4.50) -EE «nA„e,tW + (^*EE «nA„e±)(t)-0. 

± n=l ± n=l 

The operator — / + iV* is invertible, which implies 

no 

(4.51) 5^5^ a±A„e±(0 = 0. 

± n=l 

Transition from (|4.48l) to (j4.5ip shows that similarly 

no 



(4.52) «nA^e±(t)=0, fc = 0, 1, no - 1 . 

± ri=l 

The notation 

Vn^Y.^'n^tit), l<n<M 

± 

allows to view the series of conditions (|4.52p as a linear algebraic system with Van- 
dermonde matrix {Aj^}, < k < no — 1, 1 < n < uq . Hence, all rjn — . Linear 
independence of the functions and (for every index n) implies = . Yet, the 
finiteness of the sum (no < oo) is required for this conclusion. 
Step B. We prove that the equality 

(4.53) EE«»e"W=0 

± n>l 

is not possible for a nontrivial sequence {a^} G where convergence is understood 
in the £^(0, T)-norm; T > Tq . 
We introduce the function 



(4.54) ^w = -EE«« 

± n>l 

Since {a±} eP, Ee L^{Q,T). 

Step Bl. We prove that a„ = 6,i/a;„, where {6„} e P . It is assumed here that 
none of Wn = 0, and this assumption is not restrictive. For example, we may divide 
bn by n instead of w„. 

According to (|i7^ and (|T55)) . 

(4.55) £;(*)= ^^a±<(i). 

± n>l 

Using {a^} G P yields |a^|/fi < oo • Since |i?^(i)| < cq/u (see (|4.26p ) we conclude 
by the Weierstrass Af-test that E{t) £ C[0,T]. Also by Lemma 2(c), E^{0) = for 
every n, so that E{0) = . We further compute the derivative of E{t) given by (|4.55p 
term- wise and then prove that the resulting series converges in L^{0,T). We find 

(4.56) i?'(t)^EE«"(^")'W- 

± n>l 
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Here and below the symbol ~ denotes the formal differentiation of a series. We study 
this series with the help of representation (I4.36P (see Lemma 3). Substituting the term 
Qn,i{t) into the series yields a uniformly convergent series since the sequence {a„} € P. 
Substituting any of the remaining exponentials produces a series that converges in 
L^{0,T) since the set of exponentials forms an £-basis in L^(0,T), T > Tq. 
Hence, E{t) given by (|434l) or (|435l) belongs to H^{0,T), so that 

5]^a±e±-*ei/i(o,r). 

± n>l 

Since the family {e±*""*} forms an £-basis in L'^{0, T), T > Tq, the family {e=^*""7w„} 
forms an £-basis in H^{0,T) gS], [3 Sec. 11.5]. We conclude that {a^w„} G P or 
at = bt/i^n, where {bt} & P ■ 

Step B2. We prove that a„ = c„/A„, where {c„} S P (we remind here that 
A„ — uj^). Again, differentiating the series (|4.55|) twice (formally) yields 

(4.57) ^^"w-EE -(Etnt)- 

± n>l " 

We study this series with the help of representation (|4.45p (see Lemma 4) , according to 
which the structure of the ratio {Et)"{t)/ujn is the same as the structure of {Et)'{t). 
We conclude that the series ()4.36p converges in i^(0,T). 

We conclude that E G H^{0,T). Representation (|4.54p now shows that 

(4.58) E{t) = - E E 6^*'""*+"* e H\0, T) . 

± n>l 

Hence, {a^A„} £ P or at — c„/A„, where {c„} G P . Since we may now differentiate 
twice the series (I4.53P term-wise, we can repeat the proof from Step A above. We 
conclude that all at = 0, and this completes the proof of cj— independence. • 

Analyzing the proofs in Sections 4.1 and 4.2, we see that our results can be 
formulated as follows. 

Theorem 4. Let the family |e*'^"*} forms an C-hasis in L^{0,T), and functions 
en(t) satisfy the equations 

(4.59) en{t)+filen{t)+nl [ N{t ~ T)en{T)dT ^ 

Jo 

with the initial conditions 

(4.60) e„(0) = l, e„(0) = ^Ain• 
Ti2e^ the family {e„(t)} is aiso an C-hasis in L^(0, T) . 

To prove this statement we notice that the £-basis property of {e*^"*} implies 

that 

(4.61) sup5|^„|<oo, |^„| + 1 >C-n 

with some positive constant C (see, e.g. [B] Sees. H.l, H.4]; the addition of 1 in the 
last inequality reflects the possibility of = for one single n.) In turn, the proofs 
in Section 4.1 demonstrate that conditions ()4.6ip imply the quadratic closeness of the 
families {e„(t)} and {eV.t+JV(o)t/2| 

Theorem 3 follows now from Theorem 4 and Proposition 2. 
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5. Controllability. Observability. Now we apply the basis property results 
(Theorem 3) to the moment problem (j3.13p which we may rewrite after simple change 
of notations in terms of the family {e^{t)}, n > 1 as follows 

(5.1) a± / 9{t)et{t)dt 

Jo 

with {a±} e £2 _ 

Theorems 1 and 2 follow now from the results of [6^ (Theorems I.2.1(a,e), 111.3.10(a), 
and Lemma III. 2. 4). We can summarize the facts we use from the cited references as 
follows. 

Since the family {e^{t)} forms an £-basis in L^{0, T) for T > Tq, then 
(i) if T > To, the moment problem (15.1)) is solvable for any {a^} e i^; 

(a) for any T > 0, the function t !—> ^Jq g{T)e^{T)dT^ is continuous from [0,T] 

toe. 

Since for T < T^, the family £ \ £^ forms a Riesz basis in L^{0, T) , then the map 
g I— >■ {a^} from L^{0,T) to (} defined by (15.1)) has the range with codimension equal 
to the number of elements in ST- 

As usually, the controllability result can be presented as observability of the dual 
system. Using the standard techniques based on the integration by parts, one can 
check that the dual system to ()2.ip - ()2.3p can be written in the form (after the change 
of variables T — t ^ t) 

(5.2) p{x)vu{x,t) = {Av){x,t)+ f N{t-T){Av){x,T)dT, {x,t) € (0,/) x (0,r), 

(5.3) v40,t)^v{l,t)=^0, 

(5.4) v{x,0) = vo{x), vt{x,0) = vi{x); vq € "H, ui e H-i . 

The equivalent form of Theorems 1 and 2 reads as follows. 

Theorem la. Let Assumptions 1 and 2 be satisfied. Then for any T > the 
estimate 

(5.5) IKO,-)lli.(o,T) <c(IKII?i + IKIlLJ 

holds with a positive constant C independent of vq, vi. 

Theorem 2a. Let Assumptions 1 and 2 be satisfied. Then for T > Tq the 
estimate 

(5.6) \HO,-)\\Uo,T)>C{\\vo\\li + \\vi\\liJ 

holds with a positive constant C independent of vq, vi. 

Equivalence of Theorems 1 and la and, correspondingly. Theorems 2 and 2a 
expresses the standard relations between a linear operator and its adjoint in a Hilbert 
space. Specifically for partial differential equations, relations between Theorem 1 
and la are described by the transposition method |30l. The estimate ()5.6p presents 
the observability inequality for the system ()5.2)) - ()5.4p which is equivalent to exact 
controllability of ()2.ip - (l2.3p . More details about applications of the duality principle 
to hyperbolic equations with memory can be found in |31| . 
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